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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3130 


A PROCEDURE FOR THE DESIGN OF AIR-HEATED ICE-PREVENTION 

SYSTEMS 1 
By Carr B. Neel 


SUMMARY 


A procedure proposed for use in the design of air-heated systems 
for the continuous prevention of ice formation on airplane components is 
set forth. Required heat-transfer and air-pressure-loss equations are 
presented, and methods of selecting appropriate meteorological conditions 
for flight over specified geographical areas and for the calculation of 
water-drop- impingement characteristics are suggested. 

In order to facilitate the design, a simple electrical analogue was 
devised which solves the complex heat-transfer relationships existing in 
the thermal- system analysis. The analogue is described and an illustra- 
tion of its application to design is given. 


INTRODUCTION 


The first designs of aircraft ice-prevention systems utilizing heated 
air were based on the specification of an arbitrary surface-temperature 
rise in clear air for the wing and tail surfaces and of an arbitrary heat- 
flow requirement through the windshield outer surface (e.g., refs. 1, 2, 
and 3) • Subsequently, flight research conducted in natural icing condi- 
tions (refs. 4 to 8) provided information enabling a more rational approach 
to the design problem whereby a system could be designed on a wet-air 
basis. Along with the research on heat requirements, sufficient data 
also have been obtained on the meteorology of icing and the rate and area 
of impingement of water drops on airfoils so that an ice-prevention system 
can be designed which will provide reasonably adequate protection during 
flight over arbitrarily selected routes. 

This report illustrates a method for the design of ice-prevention 
equipment in which the components are continuously protected by means of 

x This report was prepared originally as one of a series of lectures pre- 
sented at the Airplane Icing Information Course sponsored by the 
University of Michigan, March 30 - April 3 , 1953 ♦ 
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Internally circulated heated air. Use is made of existing information 
obtained through previous research. In addition, apparatus is described 
and demonstrated which should aid materially in the thermodynamic design 
of an air-heated system. The heat-transfer and air-flow data presented 
in this paper are applicable only at subcritical speeds; however, in gen- 
eral, the design procedure delineated herein should be usable at all speeds. 


DESIGN PROCEDURE FOR WING AND TAIL SURFACES 

Requirements for Adequate Protection 
in Icing Conditions 


The basic requirement for wing thermal ice -prevention equipment is 
that the wing must be maintained sufficiently free of ice accretions to 
enable continued safe flight in icing conditions. Some indication of the 
amount of ice, either in the form of primary accretions or runback, a wing 
can accumulate without deleterious effects is given by the tests reported 
in reference 9 . Although under certain conditions, some tolerance for ice 
formations is shown, the recommended procedure at the present time is to 
design for complete evaporation of all impinging water in the most severe 
conditions likely to be encountered in a region of continuous icing. A 
wing designed on this basis will remain clear of ice under most situations 
and will be capable of operation in the most severe conditions likely to 
be experienced, with only a relatively small amount of runback forming. 
Furthermore, such conditions are encountered only rarely and for short 
periods of time. 


Establishment of Tentative Heated-Wing Design 


Before analysis of a thermal system can be undertaken, a tentative 
design must be established. Due to the other requirements associated with 
airplane design, there may be little latitude in the thermal-system design, 
and the analysis may be resolved into a determination of the required rate 
of heated-air flow for a set configuration. Several typical designs of 
wing systems and their advantages will be discussed. 

The most efficient and generally accepted designs utilize a double- 
skin arrangement which directs the heated air across the inner surface of 
the wing. Usually, the air is delivered to the double-skin gap through 
a spanwise duct formed by the back of the inner skin and a baffle located 
several inches behind the wing leading edge. The air enters the gap near 
the stagnation line and is directed aft through the double-skin passages. 

A number of different inner-skin arrangements have been used, severai of 
which are shown in figure 1 . 


Due to the roughness presented to the air in the spanwise duct by the 
back surfaces of the corrugated-type passages, it is often desirable to 
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install a liner against the backs of the corrugations to allow a smoother 
duct air flow (see fig. 1). In this case, a gap is left either at the 
rear of .the liner or at the leading edge to allow the air to enter the 
corrugations. When the gap is located at the rear of the liner, the air 
is forced to flow forward against the portion of the corrugation joined 
to the outer skin, then aft inside the corrugation. This arrangement pro- 
vides, in effect, a double-pass heat exchanger, which tends to concentrate 
the heat in the leading-edge region. A disadvantage of this system is 
that the flow of heated air is considerably more restricted than if the 
gap were located at the leading edge, due to the added length and tortuous - 
ness of the path. In some designs, it may be necessary, structurally, to 
retain the wing nose ribs. In this case, a liner to direct the air through 
the nose-rib openings, similar to the corrugation liner, undoubtedly will 
be required. 

In efficient heated-wing designs, the heating should be concentrated 
as much as possible in the area of water-drop impingement to obtain evap- 
oration of all intercepted water. Introducing the heated air into the 
double-skin passages at the leading-edge line of the wing, therefore, is 
the best arrangement. The use of a liner to aid in concentrating the heat 
in the leading-edge region, as noted previously, may have some advantage 
in this regard. In all cases of double-skin construction, there is a con- 
duction of heat from the inner skin to the outer at the points' of attach- 
ment. Advantage is taken of this conduction in the corrugated-type inner 
skin, where the area of attachment is large, and is minimized in the case 
of a dimpled skin, where the area of attachment is small. 

The chordwise extent of the double-skin area is determined by the 
maximum distance of water-drop impingement. Usually, the required 
percentage -chord coverage will vary from the wing root to the tip. 


Selection of Airplane Operating Conditions 


The establishment of flight conditions for the design of a thermal 
ice -prevention system must be made bearing in mind the intended use of the 
airplane. An airplane intended to be flown at high altitudes, say over 
mountainous terrain, will have different requirements from one with a rel- 
atively low service ceiling. In general, since low-altitude flight is 
required of all airplanes for approach and landing and high altitude often 
is critical for an air-heated system, several operating conditions, such 
as climb, high-altitude cruise, and letdown, should be chosen for analy- 
sis to determine the critical condition. 

In selecting an airspeed for design, consideration should be given 
to the effects of speed on the heating requirement. The rate of removal 
of heat from the outer surface is increased with increase in speed, so 
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that high speed becomes a critical condition for an air-heated system. 
Nevertheless, maximum cruising speed should not necessarily be chosen as 
the design condition. Selection of the design speed should also be guided 
by the fact that it is common practice for pilots to reduce airspeed when 

flying through the turbulent conditions which frequently accompany airplane 
icing. ‘ ' 


Selection of Design Meteorological Conditions 


Because of a lack of sufficient data, selection of meteorological 
conditions for design in the past has been somewhat arbitrary. A recent 
probability analysis of icing conditions (ref. 10), however, provides data 
from which a designer can choose a wide range of simultaneous combinations 
of meteorological variables, each having the same probability of being 
exceeded. Since three variables are involved ( liquid -water content, drop 
size, and air temperature), a three-dimensional plot is obtained which 
contains an equiprobability surface, such as that shown in figure 2. By 
use of this equiprobability surface, the most critical combinations of 
meteorological variables can be established for the thermal system. 

As was mentioned previously, conditions representative of the max imum 
to be encountered in a continuous -icing situation are considered reasonable 
design criteria. The data presented in reference 10 for layer-type clouds 
are representative of continuous icing. In selecting design data from 
reference 10, the designer must choose an exceedance probability, that is, 
the number of icing encounters which the airplane might be expected to 
experience in order to exceed the design condition once.. Likewise, the 
anticipated horizontal extent of the icing clouds must also be determined. 
Suppose, for example, that the airplane is to be operated in the plateau 
region of the United States (see fig. 3) > the system is to be designed 
for an exceedance probability of 1 in 100 icing encounters, and the hori- 
zontal extent of the icing condition is chosen as 50 miles. The equiprob- 
ability surface for this situation would appear as shown in figure 4, in 
which the surface is represented as a two-dimensional plot with constant- 
temperature contours. The thermal system, then, must provide protection 
for all combinations of liquid-water content, drop size, and air tempera- 
ture presented in this figure. 

As a means for rapid evaluation of the critical combinations of the 
meteorological variables, it has been suggested that the exceedance proba- 
bility for a given component be evaluated considering water-interception 
rates alone, rather than a n alyzing the thermal-system performance for all 
combinations of water content, drop size, and free-air temperature (refs. 10 
and 11) . Such an approximation is possible in the case of a heated wing 
because the heating requirement for evaporation of all impinging water 
depends largely on the amount of water collected. Therefore, by construct- 
ing curves of constant weight rate of water-drop impingement on the 
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Ch f a \ illustrate * ^ figure 4, the critical combina- 
' °£ iijnid-water content, drop size, and air temperature are estab- 
lished for the section of wing under consideration. These curves are 
determined from knowledge of the water-drop trajectories. Each curve of 

TZltT. A ° f /f T € 2Ptl ° n# th6 ^ defines > *<”#**> the Performance of 

^°f-. heate i, WinS designed to evaporate all impinging water at a 
given condition. The critical condition is denoted as the point at which 
therate-of-interception curves and the equiprobability curve become tan- 
gent. In the example presented in figure 4, the critical icing condition 
for an air temperature, say, of 20° F is 0.24 gram per cubic meter at 
0 microns^drop diameter. At this point, the critical rate of water 
impingement is 4.3 pounds per hour per foot span. Since the probability 
ncountering high values pf water content and, hence, high rates of 
impingement, is much greater at high air temperatures than at low tempera- 

curves for the higher air temperatures represent the more crit- 
ical icing conditions for wing systems. 


Determination of Rate, Area, and Distribution 
of Water-Drop Impingement 


^ lnG f the rate > area > and distribution of water-drop impingement on 
the wing have an important influence on the heating reqSreJSntf ?heJe 

° rs mast be evaluate d to enable proper thermal-system design. The 
° f ^ aJe ^ tor f es of water approaching various shapes in flight 

of th the objecb of a number of investigations (refs. 12 to 28). Mc£t 
of the procedures developed for the calculation of trajectories are quite 

of the°airfon req ? 1 r' ? addition to knowledge of the flow field ahead 
of the airfoil, elaborate computing equipment, such as a differential 

analyzer. As a result, the cylinder-substitution procedure (ref. 6) has 
been largely used in the past for the estimation o? impingement on air! 

fref S *29? 0W Dartl 13 ^accurate in LSces 

(ref. 29), particularly for the case of low-drag airfoils. 

Two methods have been developed recently which enable a more accurate 

need^or^^differentla!^ 61 ""*^^” 1111 ^^ 1186111611 ^ charac teristics without the 
need for a differential analyzer or for laborious calculations. The first 

of these (ref . 16) is by far the simpler procedure, requiring only a knowl- 
ege of the airfoil velocity distribution to establish the importLt 
impingement parameters. A possible limitation of the methodL that it is 
based on data obtained from fairly thick airfoils (15-percent chord) and 
there is some question as to its applicability to thL aLf oil section!! 
However, in the absence of* definite knowledge of the liTnU-nt-inno 4.-., 
cednre is suggested fox costing the 

e^ennage surfaces, provided, of course, that more exact data Lfn?t 

V +vi^! b1 / ^ thS particular case under consideration. The second of these 
methods (refs. 15 and 23) requires the construction of some equipment but 
once assembled, the apparatus enables the rapid and acc^^e elation ' 



6 


NACA TN 3130 


of trajectories for any two -dimensional shape. This system appears very 
attractive for the designer confronted with provision of icing protection 
for a variety of configurations. 


Heat-Transfer Relationships 


External heat transfer and evaporation .- During flight in icing con- 
ditions, a heated wing is cooled by convection, by evaporation of the water 
on the surface, and, in the region of water-drop interception, by the water 
the wing. The rate at which heat must be supplied in order to 
maintain the wing surface at a specified temperature is, therefore, a func- 
tion of the rates of convection, evaporation, and water impingement. The 
effect of radiant heat transfer is negligible. 

The unit heat loss from a surface exposed to icing conditions and 
heated above 32° F is expressed by the following equation (all symbols 
presented herein are defined in Appendix A) : 


q = h a (t s - t d ) + h a (X - l)K(t s - t d ) + M w c w (t s -.tv) (1) 


in which the expressions on the right-hand side represent the losses due 
to convection, evaporation, and water impingement, respectively, and where 
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The rate of evaporation, or mass transfer, from the heated surface 
is given as 
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( 2 ) 


The above equations were taken from reference 30. A simple graphical 
solution to equations (l) and ( 2 ) has been devised in reference 30 , ena- 
bling quick and accurate calculations to be performed in the analysis of 
ice-prevention equipment. Equations similar to those shown above, but 
with a somewhat different arrangement of variables, are given in refer- 
ence 31 , which likewise presents graphical solutions. 

In order to solve equations (1) and (2), knowledge of the magnitude 
of the convective heat-transfer coefficient h a is required. The most 
accurate methods for evaluating this coefficient (ref. 32 ) depend on a 
determination of the boundary-layer characteristics, a somewhat lengthy 
procedure. A simplified method sufficiently accurate for design purposes 
involves the substitution of equations for cylinders and flat plates to 
approximate the heat -transfer coefficient over a wing (refs. 33 and 3 ^) . 
These equations are as follows: 


Leading-edge region 


0.19^ T f °* 49 


(¥)°" s 


1 - 



( 3 ) 


Laminar boundary layer (beyond leading edge) 


*a 


= 0.0562 Tf 0 * 5 



(*) 


Turbulent boundary layer (beyond leading edge) 


h a = 0.51 Tf 0 * 3 dill 

gO.2 


0 .8 


( 5 ) 


In the above equations, 7 represents the specific weight of the air at 
temperature Tf and the prevailing pressure. A graphical presentation 
permitting rapid evaluation of equations ( 3 ) , ( 4 ) , and ( 5 ) is contained 
in reference 30 . 
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The most uncertain factor influencing the convective heat-transfer 
coefficient is the location of transition from laminar to turbulent flow 
in the boundary layer and the chordwise extent of the transition region. 
Disturbance of the boundary layer by the presence of water on the wing 
surface undoubtedly causes premature transition. Experimental evidence 
that the start of transition under such conditions occurs near 
the rear extremity of water-drop impingement (ref. 7). The assumption 
°f start of transition at this point appears to be a good approximation. 
At the present t ime } the extent of the transition region can only be esti- 
mated. Based on limited information (refs. 5 and 7 ), it is suggested that 
the transition region be assumed to extend chordwise along the surface for 
a distance roughly equal to three-fourths of the distance from the stag- 
nation point to the start of transition. Although the extent of the tran- 
sition zone may possibly exceed this value, its use will be conservative. 

In the transition zone, a linear variation of heat-transfer coefficient 
with distance may be taken. 

A second quantity required in the solution of equations (l) and (2) 
is the wetness fraction K. In the area of impingement, measurements 
indicate a value of 1. Aft of the impingement area the water runs back 
in rivulets, and flight and tunnel data show a rapid decrease from 1 to 
about 0.3, with a gradual decline in the value of K thereafter. The 
measured variation in wetness fraction aft of the area of impingement is 
shown in figure 5, which was taken from reference 35. 

Internal heat transfer .- In the transfer of heat from the heated air 
inside the wing, the outer skin receives heat directly from the air in the 
double-skin passages and, indirectly by conduction from the inner skin 
through the joints between the skins. The inner skin is heated by the air 
in the passages and by the air behind the passages in the leading-edge 
duct. In cases where a duct liner is installed, the inner skin is heated 
by conduction through the common joints, and the liner, in turn, receives 
heat from the air in the duct. If the liner is made of an insulating 
material, this conduction very likely will be negligible. 

The transfer of heat from the air in the double- skin passages of 
heated-wing systems is unique in that the usual equations for heat trans- 
f®!" from air flowing inside tubes do not apply. Due to severe entrance 
conditions and curvature of the passages, the air flow apparently is 
entirely turbulent, and the heat-transfer coefficient is highest at the 
gap entrance, decreasing to a constant value downstream (refs. 36 and 37). 
For typical heated-wing designs, the equations derived by Hardy and Morris 
are considered the most representative. These are 




= O.I52 


(6) 
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for the entrance region, and 


h ap d h 

k “ 

beyond the entrance region when 
The entrance length extends for 25 hydraulic diameters from the entry to 
the passages and, in this region, equation ( 6 ) applies. Beyond this point, 
the value of ^ap ds constant and can be computed by use of equation ( 7 )* 

Equations ( 6 ) and (7) were determined empirically from tests of a corru- 
gated double-skin arrangement (ref. 36). It is suggested that, when a 
duct liner is employed with the air-entrance gap located at the rear of 
the liner (fig. 1 (b)), the entrance length for equation ( 6 ) be measured 
from this point. 

Measurements of the heat conduction through the joints of double- 
skin assemblies are comparatively few. References 3 6 and 38 contain 
results of measurements of the thermal resistances of typical double-skin 
joints. Reference 38* in particular, gives thermal-resistance values for 
five different assembly procedures, encompassing riveting, cementing, and 
spot welding. Results of more general investigations of the resistances 
of joints are presented in references 39 and 40 . The measurements reported 
in references 36 and 38 are in good agreement, and the values given in 
these two reports are considered the most representative presently avail- 
able. When a leading-edge-duct liner is used, these data may also be 
applied in the calculation of heat flow from the liner to the inner skin. 

An analysis of the data of reference 36, which presents limited 
measurements of the transfer of heat from the air in the leading-edge 
duct behind corrugated double-skin passages to the inner skin in the 
absence of a liner, indicates that the following equation applies: 


0.025 (!i hV ‘ 8 


( 7 ) 


steady pipe flow has been established. 


had. d h 

k 


0.30 


Gd h \ o .6 

~) 


( 8 ) 


wherein the term d^ represents the hydraulic diameter of the leading- 
edge duct. 
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In the case of a smooth inner skin utilizing spacer strips (fig. 1(a)) 
or, in the case of a duct liner, the value of h^ may be determined from 
the equation for turbulent pipe flow, as suggested in reference 5: 



(9) 


Heated-Air Pressure Losses 


Pressure .loss in the double-skin passages .- As in the case of heat 
transfer in the double-skin passages, the usual methods for calculating 
pressure losses are not applicable. The effects of curvature and entrance 
condition apparently cause higher pressure losses than predicted. The 
data considered the most representative for heated-wing systems showed no 
difference in the pressure-drop measurements for isothermal or nonisother- 
mal flow (ref. 37). In evaluating these data to obtain a generalized 
relationship, the conventional equation involving entrance head loss and 
friction loss was employed: 



( 10 ) 


where N is the entrance coefficient and f is the friction factor. 
Figure 6 presents values of friction factor as a function of Reynolds 
number, as computed from the data of reference 37 using the suggested 
value of 0.75 for Although these data were obtained with one spe- 

cific type of double skin (corrugated), it is proposed that pressure 
losses for other types of double skin be computed from equation (10), 
using N = 0*75 and values of f from figure 6. 

It will be noted that the pressure loss in equation (10) and in sub- 
sequent equations is given in terms of the hydraulic diameter rather than 
the hydraulic radius, as used in reference kl; hence, the friction factors 
presented in this report are four times those appearing in reference kl. 

Pressure loss in the leading-edge duct .- The pressure drop of the 
air flow in the leading-edge duct may be calculated from equation (22c), 
page 128 of reference kl, or 
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G 2 Av fLG 2 v av 
~ + 2gd h 


(lla) 


where Av represents the change in specific volume of the air in the sec- 
tion of duct under consideration. In most cases, Av is negligible and 
equation (lla) reduces to 


AP 


fXjG 2 V a y 

2gdh 


(lib) 


For smooth ducts, such as would exist in designs using a smooth inner skin 
or duct liner, the value of the friction factor f is given by 


0.18k 

(Od^/n ) 0 - 2 


( 12 ) 


which is a foim of equation (9a), page 119, reference hi. When a 
corrugated-type inner skin is used without a liner, the value of f for 
the duct may be approximated from 


f> _ 0.22 

' (Gdh/d)^ 


(13) 


which is representative of the data obtained for flow in rough pipes . 
(fig. 91, P* 118, ref. 4l) . 


Calculation of Required Rate of Heated-Air Flow 


The calculation of the heated-air-flow rate through the double-skin 
passages necessary to evaporate all impinging water under the critical 
design condition requires the solution of a very complex heat balance. 

The rate of evaporation and dissipation of heat from the outer surface 
is a function of the surface temperature> which, in turn, is a function 
of the heated-air -flow rate and temperature. The temperature of the air 
at any point in the double skin depends on the amount of heat previously 
removed from the air. In addition, conduction of heat from the inner to 
the outer skin and chordwise along the skins complicates the picture' even 
further . 
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One method used previously in thermal -system analyses (refs. k2 and 

ut ^ lizes tlle simplifications of assuming a constant surface temperature 
throughout the heated area and neglecting heat conduction. The heated-air- 
flow rate required to deliver the resulting quantity of heat is then com- 
puted. This procedure has the disadvantage that errors of an unknown 
magnitude are introduced, due to the fact that the actual surface temper- 
ature rarely is uniform and heat conduction in the double-skin system can 
have an important influence on surface temperature. A more accurate method 
for calculating the required air-flow rate is to divide the double skin 
chordwise into several segments and analyze the system by use of a step- 
by-step trial-and -error procedure (ref. 3 ). A heat balance is obtained 
for each segment in sequence as the passage of the air through the double- 
skin gap is analyzed. Such a procedure is extremely laborious, especially 
in view of the fact that in a heated-wing analysis involving variable sur- 
face temperature, the effective coefficient of external heat transfer 
varies with the assumed surface temperature. In addition, both upper and 
lower surfaces of the wing at several spanwise stations must be analyzed 
to establish the required heated-air-^flow distribution. 

In order to facilitate the calculation of required flow rate, the use 
of an electrical-analogue system to simulate the flow of heated air in the 
wing can aid considerably. In utilizing the electrical analogy, an elec- 
tric model is assembled to represent a portion of a wing heated inter nall y 
with air. Such a model was constructed and used to illustrate the analysis 
of a hypothetical air -heated wing. The model, which is shown in figure 7 
consists of an arrangement of electrical resistances and capacitances con- 
nected in such a manner as to simulate a thermal circuit. In the electri- 
cal circuit, electrical resistance and capacitance represent thermal 
resistance and capacitance, current flow represents heat flow, and voltage 
difference represents temperature difference. The resistances of the ana- 
logue network consist of rheostats and fixed resistors, whereas the capaci- 
tances consist of condensers. A rotating commutator and slip-ring assembly 
is employed to switch the condensers through the circuit to simulate the 
travel of heated air through the double-skin passages. Each condenser 
receives an electrical charge immediately before it is switched into the 
circuit. This represents the initial heat energy contained in the air 
before its entry into the double-skin gap. As the condenser is switched 
through the circuit, part of its charge is lost, representing the loss of 
heat of that portion of the air and its corresponding decrease in temper- 
ature. A microammeter is used to measure the current flow at various 
locations. This provides a measure of the heat flows and temperatures. 

A simplified circuit diagram illustrating a portion of the analogue system 
is presented in figure 8 . The theory of the electrical analogy is dis- 
cussed briefly in Appendix B. 

The electrical -.analogue arrangement permits rapid determination of 
the required air-flow rates. It automatically takes into consideration 
such complex effects as air -temperature drop and heat conduction within 



MCA TN 3130 


13 


the double-skin system. A certain amount of trial-and-error procedure is 
still required, but, nevertheless, the answer is quickly established. 

Two basic steps are required in this trial-and-error procedure . These 
are (1) to obtain the surf ace- temperature distribution for an assumed 
heated-air-flow rate, and (2) to compute the total rate of evaporation of 
water for this temperature distribution. If the total rate of evaporation 
does not equal the rate of impingement, the procedure must be repeated for 
different air-flow rates until the calculated rate of evaporation is made 
to equal the rate of impingement. When the required flow rates have been 
determined for both upper and lower surfaces at several spanwise stations, 
the pressure losses must be calculated for the double-skin passages and 
the leading ^-edge duct to assure that a pressure balance has been attained. 
Very often a pressure balance will not be achieved such that an ideal flow 
distribution will result. In such cases, portions of the wing will either 
be overheated or underheated. If these conditions are intolerable, it 
may be necessary to investigate a different double-skin arrangement. In. 
any case, employing an electrical-analogue system will facilitate the 
heated -wing analysis. An example illustrating the application of the 
electric model to such an analysis is given in Appendix C. 


Discussion of the Use of an Electric Model 
in Heated-Wing Design 


The utilization of an electrical analogue in the design of air-heated 
wing systems for ice prevention has two important advantages over previous 
methods of analysis. First, it enables the .rapid determination of the 
correct heated-air-flow rate for any chosen double-skin configuration. 

Such an application of the analogue technique is demonstrated in Appen- 
dix C. Second, due to the versatility of the analogue, its use permits 
easy variation of configurations and conditions to guide the designer in 
improving the efficiency of the thermal' system. Two specific examples of 
this will be cited. 


In calculating the analogue resistances representative of the cor- 
rugated double-skin arrangement considered in the example of Appendix C, 
it was noted that a fair degree of insulation, in effect, exists between 
the inner and outer skins of a typical corrugated-passage design, even 
though the skins have been riveted together . As a result, tests were 
performed in which the thermal resistance of the joints was decreased to 
zero, a situation which might be approached in practice by seam-welding 
the skins together. It was found that, for the particular cases consid- 
ered, the flow rate of heated air could be reduced from 10.0 to 7.3 po unds 
per hour per passage at the wing-root station, and from 10.8 to 7.7 pounds 
per hour per passage at the wing-tip station, all other conditions con- 
stant. This represents decreases in the heat requirement of 27 anH 29 
percent. 
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Such a revision obviously increases the thermal efficiency of the 
double-skin system. This is exemplified in the curves of chordwise surface- 
temperature distribution shown in figure 9 for the two stations at the two 
joint conditions. It is seen that with the increased heat transfer result- 
ing from the improved thermal bond between skins, a higher surface temper- 
ature was obtained in the leading-edge region, even with lower air -flow 
rates, allowing lower surface temperatures to exist in the after portion. 
This provides a higher thermal efficiency since the heat is more concen- 
trated in the area of impingement where the surface is entirely wetted 
and, consequently, the evaporation rates are highest. 

Incidentally, the curves of figure 9 show how far from uniform the 
surface temperature of a heated wing is apt to be, and it is apparent that 
methods of analysis which assume a constant temperature may lead to error. 

A second possible means of decreasing the heat requirement was inves- 
tigated using the electric model. It becomes evident in the performance 
of an analysis of a wing thermal system that a considerable amount of heat 
is wasted in the heated region aft of the area of impingement, where the 
water runs back in rivulets, and only a small portion of the heat delivered 
is utilized for evaporation. This is illustrated in figure 10, which shows 
the distribution of evaporation rate for the wing-root station analyzed 
previously. Most of the water is evaporated in the region of impingement, 
even though this area is only three-eighths of the total heated area. It 
may be reasoned that if the wing surface were treated with a wetting agent 
which would decrease the surface tension of the water sufficiently to allow 
it to run back as a sheet, rather than in rivulets, the heat in this region 
would be utilized much more effectively. A test was conducted, therefore, 
to determine the magnitude of possible thermal saving under such a condi- 
tion. For the particular case of the wing -root station, the required 
heated-air -flow rate could be reduced from 10.0 to 6.6 pounds per hour per 
corrugation, representing a reduction of 34 percent, if the water aft of 
the impingement area were made to cover the surface completely. The dis- 
tribution of evaporation rate would then appear as the dashed line in fig- 
ure 10. Due to the increased evaporation rate now made possible in the 
area aft of the impingement limit, it is no longer necessary to concentrate 
the heat in the region of impingement, and the surface temperature for the 
entire heated area may be reduced as indicated in figure 11 . 

It should be noted that the magnitude of the reductions in heat 
requirement cited for these two examples apply specifically only to the 
cases considered. This is particularly true in the latter example, where 
the possible decrease in heat requirement is dependent upon the ratio of 
area of impingement to total heated area. It is probable, however, that 
similar reductions could be obtained for other cases of wing ice -prevention 
systems . 

A somewhat discouraging fact that becomes evident when analyzing wing 
ice-prevention equipment is that even for low values of liquid-water 
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content, very high heated-air-flow rates will be required at high speeds, 
unless improvements are made to currently employed systems. This, of 
course, is the reason why considerable interest has been shown in wing 
cyclic de-icing systems, wherein reductions in heating requirements can be 
achieved. However, by increasing the thermal efficiency through improve- 
ments in both the internal heat-transfer processes and the external evapo- 
rative processes, it appears that heating requirements can be reduced to 
the point where continuous air -heated ice -prevention systems, with their 
attendant aerodynamic advantages over cyclic systems, may be practical for 
high-speed aircraft. In this regard, the use of an electric model, as 
suggested before, can aid materially. Not only can it help the designer 
to obtain the most efficient operation from existing configurations, but 
it also can aid in guiding research directed toward increasing thermal- 
system efficiency. 

A desirable feature of the electric model is that it is relatively 
inexpensive to construct. On the basis of the experience gained through 
construction of the experimental model described herein,, it is believed 
that an analogue of this type suitable for analysis of a variety of air- 
heated components could be built for only a few thousand dollars. It 
appears quite reasonable, therefore, to construct an electric model for 
the specific purpose of designing an air-heated ice-prevention system. 


DESIGN PROCEDURE FOR WINDSHIELDS 


In this section, the criteria and equations required in the design of 
an air -heated windshield will be considered. The windshield system assumed 
for discussion will be the conventional double— panel arrangement consisting 
of two panels separated by a small gap through which the heated air circu- 
lates. 


Requirements for Adequate Protection 
in Icing Conditions 


The basic requirement for a windshield ice-prevention system is that 
the windshield be maintained sufficiently clear to allow adequate visibil- 
ity. To maintain a windshield free of ice accretions, it is necessary 
only to raise the temperature of the surface slightly above freezing, 
rather than to evaporate all impinging water, as in the case of wing ice 
prevention. The formation of runback ice accumulations from the windshield 
area is not considered sufficiently serious aerodynamic ally to warrant the 
excessive heating loads which would result by attempting to evaporate all 
intercepted water. 
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Selection of Design Meteorological Conditions 


The most severe icing conditions from the standpoint of windshield ice 
prevention are those occurring at low air temperatures (ref. 8). Therefore, 
a selection of design meteorological conditions should he guided by this 
fact.. An obviously important consideration in this selection is that only 
those conditions under which ice will form on the windshield should be 
investigated, inasmuch as it has been shown that windshield ice will not 
form when the water-drop size is small. In general, then, critical meteor- 
ological conditions will consist of large drop sizes at low air tempera- 
tures. 

As in the selection of design meteorological conditions for the wing 
system, it is suggested that the data of reference 10 be used in the deter- 
mination of the- critical icing condition for the windshield. In this 
regard, curves of constant rate of water interception superimposed on the 
equiprobability charts of reference 10 should prove useful in selecting 
the design condition. 


Determination of Rate, Area, and Distribution 
of Water-Drop Impingement 


In the absence of exact knowledge of the water-drop impingement char- 
acteristics for a particular windshield configuration, a substitution pro- 
cedure can be employed (ref. 8). This procedure involves the substitution 
of generalized shapes, for which water-drop trajectories have been calcu- 
lated, for the various windshield configurations. The shapes most appli- 
cable in this regard are ribbons and spheres, for which trajectory data 
are presented in reference 12. The particular object substituted depends 
on the type of windshield to be protected. For flat-plate and V-type 
windshields, the assumption is made that the rate of impingement on the 
windshield panel is equal to the rate of impingement on the projected 
area of the panel considered as one-half a ribbon (see fig. 12(a)). 

In the case of a windshield flush with the fuselage contour, it is 
assumed that the impingement would be the same as that on the windshield 
area projected onto a sphere having a diameter equal to the maximum diam- 
eter of the fuselage. This is illustrated in figure 12(b). The unit rate 
of water interception, then, may be approximated from the formula developed 
in reference 8. 


1 - COS 0J4 


My = 0.112 


( 1*0 
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where the terms E^ and. 0j^, the collection efficiency and. angle of impinge- 
ment for the sphere, respectively, are computed from the data of refer- 
ence 12. 


Heat-Transfer and Air -Pres sure -Loss Relationships 


External heat transfer .- The same heat -transfer processes discussed 
previously for the case of a heated wing in icing conditions apply to a 
heated windshield. Thus, the heat flow is defined in equation (1), which 
can be readily solved by use of the chart of reference 30. 

For the case of a windshield panel, the convective heat -transfer coef- 
ficient, as suggested in reference 8, may be calculated from equation (5), 
where S is measured from the leading edge of the panel for flat -plate 
and V-type windshields, and from the fuselage nose in the case of flush 
windshields . 

The value of the wetness fraction K should be taken as 1 for wind- 
shield calculations. 

Internal heat transfer .- The internal heat-transfer process for wind- 
shields is greatly simplified over that for wings in that the outer panel 
is heated only by the air in the double-panel gap and receives no heat by 
conduction, as in wing double-skin arrangements. The equation for deter- 
mining the internal heat-transfer coefficient is the turbulent-pipe-flow 
formula introduced earlier as equation (9) • This equation has been found 
to be applicable to the flow of air in rectangular ducts of large aspect 
ratio (p. 198, ref. kl) . 

Heated-air pressure loss .- The pressure loss of the air passing 
through the windshield gap is given by the previously presented equation, 
(11a). Values of the friction factor f may be computed from equa- 
tion (12) . 


Calculation of Required Rate of Heated-Air Flow 


In order to simplify the analysis, air-heated windshield systems in 
the past generally have been analyzed on the basis of an assumed uniform 
surface temperature. A more accurate solution will result by considering 
the variation of conditions for the entire surface, and by determining 
the required heated-air -flow rate accordingly. In this connection, the 
application of an electric model, such as that described for wing-system 
design, to the analysis of air-heated windshields will greatly facilitate 
the design calculations. The. same general procedure may be followed; that 
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is, to divide the double-panel system into a number of segments and to 
determine the air-flow rate required to maintain all points on the outer 
surface above freezing temperature. 


Discussion of Limitations Influencing 
Air-Heated Windshield Systems 


Heated-air systems for windshield ice prevention may not prove prac- 
tical for high-speed aircraft. Certain design considerations impose 
limitations which very likely could make excessively high air-flow rates 
required. For example, the strength of most windshield materials is 
reduced at temperatures above about 250° F, and for this reason the inlet 
temperature of the circulating air will be limited. The probability of 
large thermal stresses in the glass resulting from high air temperatures 
is another reason for limiting the air temperature. By so limiting the 
air temperature, small gap sizes or high air-flow rates will be necessary 
to obtain the required heat transfer, and the attendant air-pressure drop 
may be of such magnitude that a blower will be required. The use of thick 
windshield panels in military aircraft further aggravates the situation. 
Therefore, unless other factors make air-heating of windshields more 
attractive, • the use of electrical heating may prove to be more feasible 
for windshield ice prevention of high-speed airplanes. 


DESIGN PROCEDURE FOR RADOMES 

Requirements for Adequate Protection 
in Icing Conditions 


Protection of radomes against icing requires that the surface be 
maintained sufficiently clear to enable successful transmission of the 
radar signal. Whether or not complete evaporation of all impinging water 
is necessary will depend on the undesirability of runback ice accretions 
aft of the radome and on the tolerance of the radar system for water run- 
ning on the radome surface. Some measurements of the attenuation of the 
radar signal as a result of water on the radome surface are presented in 
reference 44 . 


Proposed Heating Arrangement 


In the following discussion, the radome will be assumed to constitute 
the forward portion of the airplane nose. Consideration will be given 
only to the thermal -system aspect of radome design although, in an actual 
case, attention must also be given to the effects of the thermal-system 
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design on the transmission of electrical energy. The heating system pro- 
posed for consideration will be a double-skin arrangement formed by the 
outer housing and an inner shell, which are separated by the circulating- 
air gap, as illustrated in figure 13. In this figure, the heated air is 
shown as entering the gap at the forward end of the radome and exhausting 
at the rear. This could be accomplished most easily by introducing the 
heated air into the radar antenna chamber, but such an arrangement may be 
impracticable due to possible temperature limitations of the radar equip- 
ment. If this is the case, reversal of the air-flow direction may be 
necessary, although this scheme probably would provide a less efficient 
ice-prevention system. 


Determination of Rate and Area of 
Water-Drop Impingement 


As in the case of windshields, probably the most expedient method 
for determining the impingement characteristics on radomes is to utilize 
a substitution procedure, in the event that more exact information is riot 
available. The method proposed is to assume that the impingement on the 
radome would be the same as that on a sphere with a radius equal to the 
nose radius of the radome. The data of reference 12 can be used to com- 
pute the rate and area of water interception on the sphere. Experimental 
results have shown the sphere-substitution procedure to provide a good 
approximation (ref. 44). Calculated impingement data for various radome 
shapes, which may be of aid in determining water-drop impingement' char- 
acteristics, are given in references 45 to 48. 


Heat-Transfer and Air -Pressure-Loss Relationships 


External heat transfer .- The external heat transfer can be calculated 
by use of equation (1). Values of the convective heat-transfer coeffi- 
cient may be computed by the method developed in reference 49 for stream- 
line bodies. Experimental data on the convective transfer of heat from 
an ellipsoid and from spheres, which should prove useful in establishing 
values of heat-transfer coefficient, are presented in references 5 0 and 
51, respectively. It is suggested that the position and extent of transi- 
tion from laminar to turbulent flow be taken as proposed previously for 
wing-system analysis. The wetness fraction K may also be taken as dis- 
cussed for wings. 

Reference 46 gives a thorough treatment of the heating requirements 
for icing protection of radomes, and should be helpful in the design of 
radome heating systems . An additional study of the external heat require- 
ments for radomes is contained in reference 52 . 
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Internal heat transfer .- The heat-transfer relationship considered 
representative of the internal-flow conditions is given in equation (9) 
which was also suggested for air-heated windshields. ’ 

Heated-a ir pressure loss .- The pressure loss of the air passing 
hrough the double-skin gap may be computed from equation (lla), with 
values of the friction factor obtained from equation ( 1 p ) , 


Calculation of Required Rate of Heated-Air Flow 


The procedure outlined for wings and windshields in the foregoing 
sections can also be utilized for radomes to determine the required air- 
flow rate for ice prevention. As before, an electric model can provide 
a rapid solution to this problem. In applying the electric model, divi- 
sion of the radome heating system for computing equivalent electrical 
values must be made in such a manner as to provide equal volumes of heated 
air between two arbitrarily selected streamlines for all segments, as 
indicated in figure 13. This is necessary to retain the air-flow analogy. 


COMMENTS ON THE DESIGN OF OTHER COMPONENTS 
OF THE THERMAL SYSTEM 


_ ' Th ® remalnder of this report will be devoted to a brief discussion 

of the design of other components of the thermal system. Detailed design 
data will not be presented inasmuch as this would entail an excessive 
discourse to treat all situations fully, and such information is ade- 
quately covered in available publications. 


Propellers 


Basically, the procedure for designing an air-heated propeller is 
the same as that for a wing. The external heat-transfer equations pre- 
sented for a wing apply equally to the propeller. Internally, the situa- 
tion is somewhat different in that the heated-air-flow direction in 
accepted designs is spanwise from the root to the tip of the propeller. 

The external and internal heat-transfer processes, including analytical 
and experimental studies to obtain improved internal heat transfer, are 
treated extensively in references 53 to 58. 

. The thermal requirements for a propeller should be based' on nkin- 
aining the heated area just above freezing temperature, rather than on ' 
evaporating all intercepted water, as- in the case of a wing. Evaporation 
oi all impinging water would impose exorbitant heating loads on the system. 
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and propeller efficiency losses resulting from runback ice accretions 
undoubtedly would be quite small (ref, 59 ). 

In the analysis and design of an air-heated propeller, the applica- 
tion of an electrical-analogue system should prove extremely effective, 
since the blade metal is comparatively thick and, with the addition of* 
internal fins, the heat-conduction effects, which are easily represented 
electrically, become very important. 


Air Inlets and Scoops 


Air inlets and scoops which are heated internally by air can be 
designed in much the same manner as heated wings. One possible exception 
to the procedure is that it may not be required that all impinging water 
be evaporated. This is a choice that should be determined on the basis 
of any adverse effects which might result from runback ice formations. 
Again, an electric model can be useful in the design. 


Air Ducts 


The primary concern with the heated-air distribution ducts is that 
they transport the air with a minimum of thermal and pressure losses. 
This requires good insulation and clean aerodynamic design. Several 
papers (refs. 60 to 64) have been written on the aerodynamic design of 
ducts. The air flow in ducts with sharp bends near the entrance region 
may become obstructed by the deposit of ice. Some information in this 
regard may be found in reference 65 . 

With the advent of cabin pressurization, a new strength-requirement 
problem has arisen in the design of air ducts. Collapsing of ducts 
located in pressurized areas can readily occur if attention has not been 
given to the reduction in strength of the duct material at elevated tem- 
peratures, nor to possible column failure of the duct resulting from 
insufficient end clearance for thermal expansion. 


Heat Exchangers 


The design of suitable heat exchangers of adequate capacity is a 
specialized field in itself. Many tests have been performed and many 
reports have been written on this, subject. A few of the more pertinent 
of these works will be cited. Reference 3^ probably contains the most 
complete information on the design of aircraft exhaust gas and air heat 
exchangers available at the present time. A good bibliography pertaining 
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to heat-exchanger design is also presented therein. An illustration of 
the design of a heat exchanger which was actually installed in an airp la ne 
thermal ice-prevention system is given in reference 66 . The effect of 
altitude on the output of an exhaust-gas-and-air heat exchanger is dis- 
cussed in reference 5 . The calculation of pressure losses in heat- 
exchanger systems is covered in references 67 and 68, and this subject is 
also considered in the analysis of reference 66 . 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Jan. 20, 195^ 
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APPENDIX A 
SYMBOLS 

area perpendicular to direction of heat flow, ft 2 
cross-sectional area of air duct or passage, ft 2 
electrical capacitance, farads or microfarads 
specific heat of air at constant pressure, Btu/lb, °F 
specific heat of water, 1 Btu/lb, °F 

hydraulic diameter of air duct or passage, — kAp ft 

perimeter’ 

electrical potential, volts 

water-drop-collection efficiency,' dimensionless 

saturated water-vapor pressure at corresponding temperature, in. Hg 

friction factor for air flow in duct or passage, dimensionless 

weight -rate of air flow per unit of cross-sectional area, 
lb/sec, ft 2 ’ 

acceleration due to gravity, 32.2' ft/sec 2 

effective coefficient of total heat transfer, Btu/hr, 'ft 2 , °F 
convective heat-transfer coefficient, Btu/hr, ft?, °F 
electrical current flow, amperes or microamperes 
mechanical equivalent of heat, 778 ft-lb/Btu 
surface-wetness fraction, dimensionless 
thermal conductivity, Btu/hr, ft 2 , °F/ft 
conduction coefficient of joints, Btu/hr, ft 2 , °F 
length or distance, ft 

latent heat of evaporation of water, Btu/lb 
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M ev rate of evaporation of water per unit of surface area, lb /hr, ft 2 
M ev T total rate of evaporation of water, lb /hr, ft span 


Miji 

M w 

m 

N 

n 

AP 

P 

R 

r 


r f 

S 

T 

Tf 

t 

V 

v 

W 

w 

X 

7 


total rate of interception of water, lb /hr, ft span 

rate of interception of water per unit of surface area, lb /hr, ft 2 

liquid-water content of ambient air, grams /meter ® 

entrance pressure-loss coefficient, dimensionless 

number of commutator segments in electric model 

air-pressure loss, lb/ft 2 

absolute static pressure, in. Hg 

rate of heat transfer per unit of surface area, Btu/hr, ft 2 
electrical resistance, ohms or megohms 
leading-edge radius of airfoil, ft 
temperature-recovery factor, dimensionless 

surface distance from stagnation point in stream direction, ft 
absolute temperature, °R 

arithmetic average of temperature of heated surface and free-stream 
air temperature, °R 

temperature, °F 

air velocity, ft/sec 

specific volume of air, ft 3 /lb 

weight, lb 

weight-rate of air flow per passage, lb /hr 
Hardy's evaporation factor, dimensionless 
specific weight, lb /ft 3 
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0 M central angle, measured from stagnation point to farthest point of 

impingement on the surface of a sphere, deg 

d viscosity of air, lh/sec, ft 

9 central angle measured from stagnation point of airfoil-leading- 

edge cylinder, deg 

w commutator rotational speed, rpm 

Subscripts 

Except where otherwise defined, the following subscripts apply: 
a external air 

a<i air in wing leading-edge duct 

ap air in double -skin passage 

av average 

d heat -transfer datum 

l local 

o undisturbed free stream 

s surface 


w 


water 
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APPENDIX B 

THEORY OF ELECTRICAL ANALOGY FOR THE SIMULATION OF HEAT 
FLOW, PRACTICAL DESIGN CONSIDERATIONS FOR AN ANALOGUE, 
AND A COMPARISON OF CALCULATED AND MEASURED RESULTS 
USING AN ELECTRIC MODEL 
Theory of Electrical Analogy 


General heat -flow analogy .- The flow of heat by conduction and con- 
vection is analogous to the flow of electricity through an electrical 
resistance. The rate of heat flow varies with the thermal resistance and 
temperature gradient, whereas electrical current flow varies with the elec- 
trical resistance and voltage difference. In both cases the variations 
are identical. Likewise, the variation of heat flow with time from a 
thermal mass is identical to the current-flow variation with time from an 
electrical capacitance. Therefore, by substitution of an electrical cir- 
cuit for its corresponding thermal circuit, heat-flow studies can be made 
by comparatively simple electrical measurements . Great flexibility is 
achieved through the ease with which values of the, electrical counterparts 
can be varied to investigate different thermal effects. The primary use 
of the electrical analogy is in the study of transient heat-flow problems 
where thermal capacitance as well as thermal resistance is of importance. 


In the electrical analogy, the following equivalent relationships 
apply: 


Thermal 


Electrical 


1 Btu 
1 Btu/sec 
1 °F 

1 °F sec/Btu 
1 Btu/®F 


1 coulomb 

1 coul./sec = 1 ampere 
1 volt 

1 volt sec/coul. = 1 ohm 
1 coul./volt = 1 farad 


In use, the values of electrical resistance and capacitance take the form 


t, 3600 L 
E '~ 


(15a) 
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for conduction, or 



(15b) 


for convection, or convection plus evaporation, and 


C = (Vol) 7Cp ( 16 ^ 

where 

A area perpendicular to direction of heat flow, ft 2 

Cp specific heat, Btu/lb, °F 

C electrical capacitance, farads 

H effective heat-transfer coefficient, Btu/hr, ft 2 , °F 

k thermal conductivity, Btu/hr, ft 2 , °F/ft 

L distance in direction of heat flow, ft 

R electrical resistance, ohms 

Vol volume, ft 3 

7 specific weight, lb/ft 3 

Heated-air-flow analogy . - In the specific application of an electrical- 
analogue system to the flow of heated air under steady-state conditions, 
the term C represents the heat capacity of the air and R, the various 
thermal resistances. For this application, the velocity of the air must 
be considered in addition to resistance and capacitance, since the amount 
of heat removed from the air in a given distance of travel is governed by 
the time required for its passage over that distance. In this case. 


RC 


v ap 

ST 


= constant 


(17a) 


The above relationship controls the selection of values such that the con- 
densers in the analogue system will lose their charge proportional to the 
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heat loss of the air as it progresses through the passage. Equation (17a) 
is derived in the following manner. In a short length of air passage, the 
thermal energy loss of the air can be equated to the heat transferred to 
the adjacent surface, or 


wcpAtg^ - hAsCtg^ - t s ) 


where At Sp represents the temperature drop of the air. This equation 
may be rearranged to give 


wc p _ 'tap " t s 
hA-s At S p 


Both sides of this equation equal a dimensionless constant. Therefore, 
substitution of equivalent electrical terms for the thermal terms may 
be made without altering the constant. Considering the left side of the 
equation, if the length of air passage is AL feet containing (Vol) cubic 


feet, and the air velocity is 


w can be expressed as 


3600 


(Vol) 7 V- 
AL ' 


and 


= 3600 

hA s 


(Vol)7CpV ap 
~ hAgAL 


= constant 


Then, since h = H in the case of convection, substitution of equations 
(15b) and (16) into the above relationship yields equation (17a). 

When an electric model is used employing a commutator arrangement, 
such as that used for the tests described herein (see fig. 7), the com- 

Vo 

mutator rotational speed may be substituted for the term in equa- 

tion (17a), or 


RCu = constant 


where 


(17b) 


■ 60 V 

W = 




(18a) 


nAL 
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In the above equation n represents the number of commutator segments and 
AL the distance of travel of heated air during the passage of one segment 
past a brush. The value of Vap may be expressed as 


V 


_ w ay 
a P “ 36 ooa p 


and equation (l8a) becomes 


= w av 
60 ApnAL 


(18b) 


The value of C represents the thermal capacity of the volume of air 
present in a segment of air passage AL feet long. If AW represents 
the weight of air in the segment, then 


C = AWcp farads 


or 


C = AWCp X 10 s 


microfarads 


Since 


AW = 


ALAp 

v a v 


then 


ALApCp 

C X 10 6 microfarads 

v av 


( 19 ) 


If equations (l8b) and (19) are combined. 
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Cto 


CpW 

60 n 


x 10 s 


( 20 ) 


where C is in microfarads. 

When the above equations are applied to the calculation of values for 
an electric model, values of R are computed using equations (15a) and 
(15b), and values of Cto are determined from equation (20). The values 
of R, C, and to may then be adjusted for convenience, provided equa- 
tion (17b) is satisfied. If the values of R or C are to be adjusted, 
all values of R or C in the network must be changed proportionately. 

For instance, if it is desired to increase a resistance value by a factor 
of 10, all resistances in the circuit must be increased by the same factor. 

Conversion of measurements .- Temperature differences are directly 
proportional to voltage differences in the analogue system, and heat flow 
is likewise directly proportional to current flow. Since temperature dif- 
ference and heat flow are directly related by the equations 


and 


q = HAt 


q. = 



or electrically, 



then the value of temperature difference may be obtained indirectly by 
the measurement of heat flow. The converse also is true. In the elec- 
trical analogy, then, a measure of temperature difference is obtained 
from measurements of current flow through known resistances. To convert 
the electrical measurements directly in terms of thermal units, a known 
electrical value must be assigned a specific thermal value. For example, 
a potential difference of 25 volts may be considered to represent a tem- 
perature difference of 100° F. Then an actual measurement of 10 volts is 
equivalent to a k) 0 F temperature difference. 

In the case of the application of the electric model, the voltage 
of the battery which precharges each condenser prior to its entry into the 
analogue circuit may be considered to represent the temperature rise above 
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ambient of the heated air in the leading-edge duct at the entrance to the 
double-skin passages. The surface-temperature rise, then, is obtained by 
determining the voltage across the external-air resistance. The equation 
expressing this relationship is 


At s Eg 

At ad E a<i 

where the subscript s denotes the surface condition and a^ pertains 
to the initial heated-air condition at the passage entrance. Since cur- 
rent, rather than voltage, measurements are made with the electric model, 
this equation becomes 


(21 > 

where R a is the external-air resistance. The value of R a „ is chosen 

a d 

arbitrarily, and I a ^ is the measured current flow through R a -, when the 
battery voltage Ea<i is impressed across this resistance. 


Practical Considerations Regarding the Design 
and Use of an Electrical Analogue 


A few of the more pertinent aspects concerning the design and use of 
an electrical analogue are worthy of note at this point. A more detailed 
discussion of the requirements associated with the design of electrical- 
analogue systems is contained in reference 69. Further information is 
also given in reference 70, which presents a comprehensive consideration 
of the design of analogues for steady-state applications requiring only 
electrical-resistance elements . 

Requirements for condensers .- The leakage resistance of the conden- 
sers used to simulate thermal capacity in the electrical network should 
be at least of the order of 100 times the effective resistance of the 
circuit in which they are connected. Since analogue-circuit resistances 
generally are several megohms, this means that the leakage resistance 
must be very high, and for most applications electrolytic condensers will 
be unsatisfactory. Most commercially available paper, mica, and oil-filled 
condensers, on the other hand, have been found to possess suitable leakage- 
resistance characteristics. Resistance values above 1000 megohms can be 
expected for condensers of these types with nominal capacitance ratings 
up to 10 microfarads. 
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The use of a conventional capacitance bridge for the measurement of 
condenser capacities prior to their connection into the analogue circuit 
has been found unsatisfactory in most cases. The conventional bridge 
measures the capacitance in an alternating-current circuit, whereas the 
condensers are used in a direct-current circuit. Most paper and oil-filled 
condensers exhibit a characteristic when charged termed "dielectric absorp- 
tion." Instead of merely a surface charge forming on the condenser plates, 
the dielectric apparently becomes permeated with the charge also. There 
is some resistance to the permeation process which increases the time for 
charging and discharging of the condenser. The effective dir ect-current 
capacitance, therefore, must be determined from charging or discharging 
curves for each condenser when connected in an RC circuit. The curves° 
should be checked against theory inasmuch as a large degree of dielectric 
absorption can render a condenser useless for analogue applications. 

Requirements for mounting panels . - The primary electrical requirement 
for the material of the panels used for mounting the electrical components 
is that it must have very high leakage resistance under all anticipated 
conditions of humidity and temperature. Tests have shown that clear 
methyl methacrylate resin is a very satisfactory material from this stand- 
point . 

Requirements for measuring equipment .- The representative temperature 
difference and heat flow at any point in an analogue circuit may be deter- 
mined from either a voltage measurement across a known resistance or a 
measurement of current flow through the resistance. There are certain 
requirements for each measuring technique. When voltage is measured, the 
resistance of the meter must be large relative to that of the analogue 
circuit. Because the resistance of most meters is comparatively low, the 
use of a cathode -follower circuit into which the meter is connected should 
provide the necessary measuring resistance. 

When current is measured, the resistance of the meter must be small 
relative to that of the analogue circuit. Generally, the circuit resist- 
ances will be sufficiently high that standard meters can be used without 
difficulty. 


Comparison of Calculated Values With Results 
From the Electric Model 


In order to obtain a check on the accuracy of the electric model for 
the simulation of the flow of heated air in a wing, a problem was placed 
on the model for which the solution had been previously computed. The 
problem considered was the flow of heated air through a circular tube 
0 inches long which was cooled externally. Surface temperatures along 
the length of the tube, as obtained by the two methods, provided the basis 
for establishing the accuracy of the analogue. 
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The heated air was assumed to enter the tube at a temperature of 
300° F, and the coefficient of internal heat transfer was taken as 
15 Btu/hr, ft 2 , °F. An air-flow rate of 10 pounds per hour was assumed. 

The tube was considered to be 0.35-inch diameter with an infinitesimally 
thin wall. Heat conduction lengthwise, therefore, was neglected. The 
external coolant was considered to be at a temperature of 0° F, and the 
external coefficient of heat transfer was taken as 40 Btu/hr, ft 2 , °F . 

In computing values of tube -surface temperature for comparison with 
the solution by the electric model, a step-by-step procedure was employed. 
The tube was divided into 1-inch-lengthwise increments, and a balance of 
the heat removed from the air to that delivered to the surface was obtained 
for each increment. 

For the electrical analogy, the tube was divided into 1-inch incre- 
ments also. Three separate tests were made with different values of 
resistance, capacitance, and commutator rotational speed. These were as 
follows : 


Test 

R int> 

megohms 

^ext> 

megohms 

c, 

mfd 

0 ), 

rpm 

1 

2.34 

0.874 

0.508 

90 

2 

1.17 

0.437 

0.508 

180 

3 

1.17 

0.437 

1.01 

90 


The three tests should give identical results, provided the above values 
are maintained rigidly, since they satisfy equation (17b) and, hence, 
represent identically analogous conditions. Results of the three tests 
are presented in figure 14, which compares the tube -surface -temperature 
distribution as measured using the electric model with the calculated 
variation. The values for the three test arrangements agree with each 
other and with the calculated curve within 1° F • From these data, it 
may be concluded that the electric model provides sufficiently accurate 
results for thermal-system-design applications. 
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APPENDIX C 

APPLICATION OF AN ELECTRIC MODEL TO THE ANALYSIS 
OF AN AIR -HEATED WIND FOR ICE PREVENTION 
Conditions and Assumptions 


Wing-geometry and internal conditions .- In order to ill ustrate the 
application of the electrical -analogue system to the analysis of an air- 
heated wing, a typical double-skin arrangement consisting of corrugated- 
type passages was chosen. The wing geometry was taken as follows: 

Root section .... NACA 6^1 -312 13-ft chord 

Tip section .... NACA 64-^08 6-ft chord 

Angle of attack .... 0° 

The analysis was performed only for the upper surface of the wing at 
the root and tip stations. Cross sections of the wing leading-edge region 
at the two stations are presented in figure 15. Details of the double - 
skin corrugated passages also are shown in this figure. The required 
chordwise extent of double-skin area was taken as 16 inches for the root 
section and 9 inches for the tip. These correspond to 8.75- and 11. 7- 
percent chord, respectively, for the two stations. Normally, the required 
chordwise extent of double skin would be calculated from knowledge of the 
largest size of water drops expected, but for this example the values 
were chosen arbitrarily. 

For purposes of analysis, the wing was divided chordwise into 8 seg- 
ments at the root and 6 at the tip, as illustrated in figure 15. Condi- 
tions were taken as being uniform throughout each segment. The analysis 
was performed for a 1-inch-wide spanwise strip encompassing a single cor- 
rugation at each wing station. 

The heated air was assumed to enter the wing-root corrugation at a 
temperature of i+OO 0 F and an absolute static pressure of 1073 lb/ft 2 . 
Assuming full dynamic pressure of the air stream is available at the air 
•inlet to the thermal system, this represents a pressure drop of 63 lb/ft 2 , 
or about 4o percent__of the dynamic pressure, up to the entrance of the 
corrugation, for the airspeed and altitude conditions given below. At the 
entrance to the wing-tip corrugation, the air temperature was taken as 
380° F and the absolute static pressure as I060 lb/ft 2 . 

In order to simplify the example, the transfer of heat from the air 
in the leading-edge duct to the back of the corrugated skin was considered 
negligible. This might be the case when a duct liner of insulating mate- 
rial is used. 
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Flight, meteorological, and water -drop- impingement conditions .- The 
assumed airplane flight conditions are 

True airspeed 350 mph 

Pressure altitude 20,000 ft 

It was assumed that the airplane was to "be operated over the plateau area 
of the United States (fig. 3). Since high values of liquid-water content 
are critical for a wing thermal system, and since high. water contents 
usually are associated with high air temperature, a temperature of 20° F 
was selected for the design condition. An exceedance probability of 1 in 
100 icing encounters was taken as the design criterion, and the values of 
water content presented in reference 10 for this condition were modified 
for a 50-mile horizontal extent. The resulting relationship of liquid- 
water content and drop size is shown in figure 4. 

In order to determine the critical meteorological condition for the 
wing, curves of constant rate of water-drop impingement for the wing-root 
section were plotted in figure 4. These were estimated from data presented 
in reference 11. In the interest of expedience, estimated curves were used 
rather than values computed from reference 16 inasmuch as these curves are 
satisfactory for illustrative purposes. From figure 4, it is seen that 
the critical condition for this section of the wing occurs at 0.24 gram 
per cubic meter and 20 microns. For purposes of illustration, it was 
assumed that this was also the critical condition for the wing-tip sec- 
tion, although it is unlikely that the same meteorological condition 
actually would be critical for both sections. The calculated total rates 
of water interception are 4.3 and 3*3 pounds per hour per foot span for 
the root and tip sections, respectively. The rearmost point of impinge- 
ment for these conditions was taken as 6 inches from the stagnation point 2 
for the root section and 4-l/2 inches for the tip section. The assumed 
distributions of impingement for the upper surface at the two stations 
are shown in figure 1 6 . In constructing these curves, all the water 
striking the wing was considered to be divided equally between the upper 
and lower surfaces. The information contained in reference 1 6 was used 
as a guide in establishing the shape of the distribution curves. 

Electrical circuit representing the air-heated wing .- The electrical 
circuit representing the thermal circuit of the air-heated wing considered 
in the application of the electric model is shown in figure 17. It will 
be noted that there are two locations of possible current measurement for 
each wing segment: one in the external-air-resistance circuit to provide 

a measure of surface -temperature rise, and one in the internal-air- 
resistance circuit to provide a measure of heated -air -temperature rise. 


2 The stagnation point and wing leading edge in this case are coincident. 
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Calculation of Analogue Values 


In determining the required values of resistance for the analogue 
circuit, it was decided to multiply all calculated values "by a factor of 
80. The reason for this is that it -was considered desirable to reduce 
the computed values of C(o by this factor to obtain reasonable values 
of capacitance and commutator rotational speed. Since the result ing 
resistances are quite large, it is desirable to convert them in terms of 
megohms. Thus, from equations (15a) and (15b) 


or 


= 0.288 L 
kA 


megohms 


0.288 

HA 


megohms 


(22a) 


(22b) 


For evaluation of capacitance and rotational speed, the specific 
heat of the heated air was assumed constant at a value of 0.243 Btu/lb, 
°F, corresponding to a temperature of 300 F. The number of commutator 
segments on the electric model is 12. If these values are substituted in 
equation (20) and a factor of l/8o is included to satisfy equation (17b) , 
the product Cw becomes 


Cw = 4.22 w 


The value of C was chosen as O.5O8 mfd, because of the availability of 
condensers of this capacity. 3 As a result, the commutator rotational 
speed is given by 


<j = 8.31 w ( 23 ) 

In establishing the required external- and internal-air resistances, 
charts depicting these resistance values were found useful. These, 
together with other charts and calculation procedures, are discussed below. 
In order to clarify the derivation of these charts, specific examples will 

Q l ' ' ' " r - ' r ' .... 

Since there are 12 commutator segments, 12 condensers were required. 
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be shown. In these examples, segment 8 of the wing-root station will be 
chosen for illustration. 

External-air-resistance chart . - The resistance to extern a l heat flow 
from the wing surface is a function primarily of the convective heat- 
transfer coefficient and the rate of evaporation of water. Because the 
rate of evaporation is a function of surface temperature, the effective 
coefficient of heat transfer varies with surface temperature, and the 
analogue resistance representing this effective coefficient must vary 
likewise. To determine this variation, curves of unit heat flow as a 
function of surface temperature were constructed for each chordwise wing 
segment using equation (1) and taking average values for conditions pre- 
vailing over the segment. Solutions to equation (1) were obtained with 
the aid of the chart included in reference 30. Values of the convective 
heat-transfer coefficient h a in equation (1) were determined by use of 
equations (3)> (Mj snd (5)« The curve of h a calculated for the wing- 
root station is presented in figure 18. Transition was assumed to start 
at the end of the area of impingement, with a gradual transitional region 
following. The value for the Wetness fraction K in equation (1) was 
taken as 1 in the impingement area. Aft of this area, the values of K 
given in figure 5 were used. 

A curve of the heat-flow variation for wing segment 8 is shown in 
figure 19 and was constructed as follows. For this segment, the following 
conditions apply: 

My = 0 lb/hr, ft 2 
h a = 36 Btu/hr, ft 2 , °F 
to = 20° F 
V 0 = 51^ ft/sec 
p, = 12.75 in. Hg 
p Q = 13.75 in. Hg 
K = 0.15 

the above values and the chart of reference 30, the heat-transfer 
q was computed for the expected range of surface temperatures. 

The effective heat-transfer coefficients at several temperatures 
obtained from the heat-flow curves by use of the relation 


From 

rate 

were 


H = 


A 


t s " id 


in which t^ represents the surface temperature when q is zero. Then, 
by substitution of these values in equation (22b), curves of R a vs. t s 
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for each segment were obtained. The curve of resistance as a function 
of surface temperature for segment 8 is given in figure 20. 

Internal-air-resistance chart .- The resistance to heat flow from the 
internal heated air to the walls of the air passages is a function of the 
weight flow of air and the distance from the passage entrance . Values of 
the coefficient hap were computed for each segment for a range of air- 
flow rates from equations (6) and (7), and the electrical resistance as 
a function of flow rate was obtained for each segment by use of equa- 
tion (22b) where H = h^. These data were computed for both the outer- 

skin surface and the inner-skin surface. 

For the case of segment 8 , equation (7) was applied, inasmuch as this 
segment is beyond the entrance region of 25 hydraulic diameters. Figure 21 
presents the internal-air-resistance curves for segment 8. In the calcula- 
tion of these curves, the values of k and d (conductivity and viscosity) 
for air were taken as constant and were evaluated at a t emp erature of 
300° F . 

Calculation of other resistance values .- The electrical resistances 
representing the thermal resistances in the aluminum skins, both inner 
and outer, were computed by use of equation (22a), where L represents 
the effective distance along the skin in the direction of heat flow. The 
value of k in equation (22a) was taken as 67.7 Btu/hr, ft 2 , °F/ft, rep- 
resentative of 2k ST aluminum alloy. 

To establish the resistance between the inner and outer skins at the 
joints, equation (22b) was used, in which kj, the conduction coefficient, 
was substituted for H. The value of kj was taken as 55 Btu/hr, ft 2 , °F, 
as suggested in reference 36 . 

Rate-of -evaporation chart ♦ - The determination of the rate of evapo- 
ration from the wing surface was facilitated by the use of curves of M ev 
as a function of t s , which were constructed for each segment with the 
aid of equation (2) and reference 30. The values of h a and K were the 
same as used in the calculations for the external-air-resistance chart. 

A plot of the rate of evaporation as a function of surface temperature 
for segment 8 is illustrated in figure 22. 

The total rate of evaporation from the entire heated area is given 
as the summation of the individual rates for each segment, or 


MeVj = Mev^L + M ev ^AL + . . . + M ev ^AL 
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Thus, 

Meyrp = AL(M eVl + M eVs + ... + M eVQ ) (24) 

The term AL represents the length of surface distance for each segment, 
and the subscript numbers designate the appropriate segments. 

Heated-air pressure-loss chart .- Curves of the parameter AP/v av as 
a function of air-flow rate for the wing-root and -tip upper-surface cor- 
rugations were computed from equation (10), and are shown in figure 23. 

The term AP represents the total pressure drop in inches of water for 
each entire corrugation. 


Procedure for Determining Required 
Heated-Air -Flow Rate 


The procedure for determining the required heated-air -flow rate util- 
izing the electric model will now be considered. To satisfy the basic 
requirement for adequate ice protection, as previously set forth herein, 
it is necessary to evaporate all impinging water under the design condi- 
tion. This is accomplished by a trial-and-error procedure involving two 
basic steps, as mentioned previously. These are (1) to obtain the 
surface-temperature distribution for an assumed air-flow rate, and (2) to 
compute the total rate of evaporation for this temperature distribution. 
This procedure will be demonstrated for the wing-root station. 

Initial estimates .- In order to set the appropriate values of 
external-air resistance on the model, it is first necessary to estimate 
the values of surface temperature, in view of the fact that the value of 
resistance is dependent upon temperature. An initial estimate of the 
approximate magnitude of temperature in the leading-edge region can be 
obtained by calculating the surface temperature required to maintain a 
dry surface. In this case. 


MvPjCp 

0.622h a 


(25a) 


where e s represents the water-vapor pressure at temperature t s . The 
contribution of ej is small and, for estimation purposes, can be neg- 
lected. This leads to 
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e s = 


For segment 1 of the root section, 

M w = 8.25 lb /hr, ft 2 
Pj = 14.8 in. Hg 
Cp = 0.24-5 Btu/lb, °F 
h a = 20 Btu/hr, ft 2 , °F 

The resulting value of e s is found to be 2.4o in. Hg, and the correspond- 
ing value of t s is 107° F. It is not necessary, however, to maintain a 
dry surface (i.e., evaporate all intercepted water in the area of impinge- 
ment) because, some evaporation will be obtained in the heated region aft 
of the impingement area. Therefore, the surface temperature of segment 1 
can be somewhat less than 107° F. As an initial estimate, a value of 
85° F was chosen. The surface temperatures for the remainder of the seg- 
ments were then estimated, bearing in mind the magnitude of the external- 
air resistances and the probable decrease in temperature of the heated 
air as it proceeds through the corrugation. A rough check of the appro- 
priateness of these values of surface temperature was obtained by computa- 
tion of the total rate of evaporation with the help of the evaporation 
charts and equation (24) . The initial estimate of the surface tempera- 
tures, together with the corresponding resistances and evaporation rates 
are tabulated below. 


Segment 

number 

Op 

R a> 

megs 

Mev 

lb/hr, ft 2 

1 

85 

0.12 

3.9 

2 

90 

.18 

2.9 

3 

95 

.22 

2.7 

4 

90 

.37 

1.3 

5 

70 

.28 

1.0 

6 

60 

.29 

.6 

7 

55 ! 

• 35 

.4 

8 

50 

.43 

• 3 


The sum of the M ev values is 13*1 and, since AL equals 2 inches, or 
0.167 foot, the total rate of evaporation, from equation (24), equals 
2.18 lb/hr, ft. Since the required rate of evaporation must equal the 
rate of impingement, or l/2 x 4.3 = 2.15 lb/hr, ft, the above estimated 
surface temperatures appear to be of the right order of magnitude, and 
the corresponding external-air resistances were used in the initial test. 


MwPZQp 

0.622h a 


(25b) 


the following average values apply: 
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The required heated-air -flow rate must he determined by a trial-and- 
error procedure in which a value of w is assumed, the resulting surface- 
temperature distribution is obtained by use of the electric model, and the 
corresponding evaporation rate is computed. The procedure is repeated, 
if necessary, until the rate of evaporation is made to equal the rate of 
impingement. 

Calculated values .- The value of heated-air-flow rate taken for the 
initial test was 10 lb/hr per corrugation. The electric model was adjusted 
for the corresponding internal -air-resistance values. From equation (23) , 
the commutator rotational speed representative of the flow rate is found 
to be 83 rpm. The initial test performed using these values produced the 
surface temperatures listed in the following table: 


Segment 

number 


First 

estimate 

Measured, 
first test 

Second 

estimate 

Measured, 
second test 

°F 

megs 

R ap> 
outer , 
megs 

R a p > 

inner, 

megs 

Atg, 

°F 

tg> 

°F 

4 S > 

°F 

R a> 

megs 

At S ; 

Op 

ts> 

°F 

M-evf 

lb /hr, ft 2 

1 

85 

0.12 

0.69 

0.58 

60 

80 

82 

0.125 

61.5 

81.5 

3.5 

2 

90 

.18 

.97 

.82 

61.5 

81.5 

85 

.20 

64.5 

84.5 

2.45 

3 

95 

.22 

l.l4 

.95 

6l 

81 

87 

• 25 

66 

86 

2.05 

4 

90 

.37 

1.17 

1.00 

78 

98 

95 

• 35 

74 

94 

1.45 

5 

70 

.28 

1.17 

1.00 

58 

78 

75 

.27 

55 

75 

1.15 

6 

60 

.29 

1.17 

1.00 

54.5 

74.5 

69 

.26 

49.5 

69.5 

• 95 

7 

55 

.35 

1.17 

1.00 

54 

74 ■ 

68 

• 29 

48 

68 

• 75 

8 

5° 

.43 

1.17 

1.00 

52.5 

72.5 

66 

•39 

46 

66 

.6 


It is seen that the measured temperatures differ from the initially esti- 
mated temperatures. Because of this, the values of resistance used in the 
test are not appropriate for the measured temperatures. The correct tem- 
peratures, therefore, lie between the estimated and measured values, and 
it was necessary to make a second estimate of surface temperature and, 
by use of the corresponding resistance values, to perform a second series 
of measurements. These latter estimates and measurements are given in the 
above table. Good agreement now exists between the measured and estimated 
temperatures, signifying that the correct values of resistance were 
employed and, hence, the measured temperatures are correct. 

For these measurements, the value of Rg^ in equation (21) was 
chosen as 0.7 megohm, and the measured value of Ig^ was 28.7 micro- 
amperes when the battery voltage was impressed across Rg^. During the 
first test, the measured value of I for segment 1 was 26.5 microamperes 
and, as noted in the table, R a was 0.12 megohm. Then, since tg^ was 
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assumed to be *400° F, Ata^ equals 380° F, and At s from equation (21) is 
6o° F; thus, t s equals 80° F. The remainder of the surface-temperature 
values were computed from the measurements in the same manner. 

The rates of evaporation for each segment were obtained from the 
rate -of -evaporation charts for the corresponding surface temperatures and 
are presented in the above table. Substitution of these values in equa- 
tion (24) yields M eVT = 2.15 lb/hr, ft, which also is the value of the 

rate of water impingement. Thus, the required heated-air-flow rate for 
the upper corrugation at the wing-root station is 10.0 lb/hr at an inlet 
temperature of 400° F. 

The measured air temperature at the corrugation outlet was I85 0 F. 

By use of this value, together with the assumed inlet air temperature and 
average pressure, the average specific volume of the air in the corruga- 
tion was computed as 39*1 ft 3 /lb. The air-pressure loss through the cor- 
rugation was determined from the pres sure -loss chart for this value of 
v av and was found to be 19.0 inches of water. 

By the same procedure, the required air-flow rate for the tip section 
was established as 10.8 lb/hr per corrugation at an inlet temperature of 
380° F. The resulting pressure drop was calculated to be 13.2 inches of 
water, which means that the pressure loss in the leading-edge duct from 
the wing root to the tip could be about 6 inches of water. 
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Figure 3* - Map of the United States showing boundaries of areas used 
in the geographical classification of icing data presented in 
reference 10. 
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Figure k.- Equiprobability surface and curves of constant rate of 
water interception for rapid evaluation of critical meteoro- 
logical conditions. 
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SURFACE -WETNESS 
FRACTION, K 



Figure 5.- Variation of wing-surf ace-wetness fraction aft of area of 
water impingement with distance from limit of impingement. 
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Figure 6.- Friction factor as a function of Reynolds number for 
determining air-pressure losses in double-skin passages using 
equation (10). 
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(a) PHOTOGRAPH OF COMPLETE MODEL 
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(b) SKETCH SHOWING DETAILS OF COMMUTATOR AND SLIP-RING ASSEMBLY 


Figure 7.- Electric model used for simulating air-heated ice -prevention 

systems. 
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22-VOLT BATTERY 
PROVIDING INITIAL CHARGE TO 
CONDENSER, REPRESENTING 
INITIAL HEAT ENERGY OF AIR 
ENTERING DOUBLE-SKIN 
PASSAGE 
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Figure 8.- Simplified circuit diagram illustrating a portion of analogue 
system for simulating the flow of heated air through a wing. 
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DISTANCE. FROM LEADING EDGE, INCHES 
(a) WING -ROOT STATION 



(b) WING-TIP STATION 


Figure 9»~ Comparison of wing-surf ace-temperature distribution for two 
conditions of conduction between inner and outer skins. 
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RATE OF EVAPORATION, 
M ev , Ib/hr, ft 2 



Figure 10.- Comparison of distribution of evaporation rate over wing sur- 
face for two conditions of surface wetness aft of impingement region. 



Figure 11.- Comparison of wing-surf ace-temperature distribution for two 
conditions of surface wetness aft of impingement region. 
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Figure 12.- Illustration of substitution procedure for calculating rate 
of water impingement on different types of windshields, using data 
for ribbons and spheres. 
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Figure 13.- Proposed heating system for radome ice prevention and division 

of radome for analysis. 
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Figure 1^.- Comparison of tube-surface-temperature distribution as cal- 
culated, and as measured using an electric model to simulate heated- 
air flow. 
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DETAIL 

(a) WING-ROOT STATION 
(NACA 64,-312 SECTION, 13-FT CHORD) 



(b) WING-TIP STATION 
(NACA 64-408 SECTION, 6-FT CHORD) 


Figure 15. - Cross sections of forward region of wing upper surface chosen 
for analysis, showing chordwise divisions for use with electric model, 
and other details. 
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DISTANCE FROM LEADING EDGE, INCHES 
(a) WING- ROOT STATION 



DISTANCE FROM LEADING EDGE, INCHES 
fb) WING -TIP STATION 


Figure 16 .- Assumed distribution of rate of water impingement on wing 

upper surface. 
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Figure 17 •- Electrical circuit representing thermal circuit of an air- 
heated wing as used in application of electric model. 



FROM LEADING EDGE, INCHES 


Figure 18.- Calculated variation of external convective heat-transfer 
coefficient along wing upper surface at root station. 
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HEAT-TRANSFER 
RATE , q , 
Btu/hr.ft 2 



Figure 19.- Heat-transfer rate as a function of surface temperature for 
segment 8 of wing surface at root station. 



Figure 20.- External-air resistance as a function of surface temperature 
for segment 8 of wing surface at root station. 
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AIR -FLOW RATE, w, lb /hr PER PASSAGE 


Figure 21.- Internal-air resistance as a function of heated-air -flow rate 
for segment 8 of wing surface at root station. 



Figure 22.- Evaporation rate as a function of surface temperature for 
segment 8 of wing surface at root station. 




